Mon. Not. R. Astron. Soc. 000, [Tl[8] (2006) Printed 1 February 2008 (MN WT^K. style file v2.2) 



Timing evidence in determining the accretion state of the 
Seyfert galaxy NGC 3783 

D. P. Summons 1 ' 3 *, P. Arevalo 1 , I. M. M c Hardy\ P. Uttley 2 and A. Bhaskar 3 

1 School of Physics and Astronomy, University of Southampton, Southampton S017 1BJ, UK 

2 Astronomical Institute 'Anton Pannekoek', University of Amsterdam, Kruislaan 403, 1098 SJ, Amsterdam, the Netherlands 
3 School of Engineering Science, University of Southampton, Southampton SO 17 IB J, UK 



Received /Accepted 



ABSTRACT 

Previous observations with the Rossi X-ray Timing Explorer (RXTE) have suggested 
that the power spectral density (PSD) of NGC 3783 flattens to a slope near zero at low 
frequencies, in a similar manner to that of Galactic black hole X-ray binary systems 
(GBHs) in the 'hard' state. The low radio flux emitted by this object, however, is 
inconsistent with a hard state interpretation. The accretion rate of NGC 3783 (~ 7% 
of the Eddington rate) is similar to that of other AGN with 'soft' state PSDs and 
higher than that at which the GBH Cyg X-l, with which AGN are often compared, 
changes between 'hard' and 'soft' states (~ 2% of the Eddington rate). If NGC 3783 
really does have a 'hard' state PSD, it would be quite unusual and would indicate 
that AGN and GBHs are not quite as similar as we currently believe. Here we present 
an improved X-ray PSD of NGC 3783, spanning from ~ 10 -8 to ~ 1CP 3 Hz, based 
on considerably extended (5.5 years) RXTE observations combined with two orbits 
of continuous observation by XMM-Newton. We show that this PSD is, in fact, well 
fitted by a 'soft' state model which has only one break, at high frequencies. Although 
a 'hard' state model can also fit the data, the improvement in fit by adding a second 
break at low frequency is not significant. Thus NGC 3783 is not unusual. These results 
leave Arakelian 564 as the only AGN which shows a second break at low frequencies, 
although in that case the very high accretion rate implies a 'very high', rather than 
'hard' state PSD. The break frequency found in NGC 3783 is consistent with the 
expectation based on comparisons with other AGN and GBHs, given its black hole 
mass and accretion rate. 

Key words: galaxies: active - galaxies: Seyfert - galaxies: NGC 3783 - X rays: 
galaxies 



1 INTRODUCTION 

Super-massive black holes in active galactic nuclei (AGN) 
and Galactic stellar-mass black hole X-ray binary sys- 
tems (GBHs) both display aperiodic X-ray variability 
which may be quantified by calculating the power spec- 
tral densities (PSDs) of the X-ray light curves. The PSDs 
can typically be represented by red-noise type power 
laws (i.e. P(v), the power at frequency v, oc u a where 
a ~ -1) with a bend or break (to a ^ -2) at a 
characteristic PSD frequency. The time-scale, correspond- 
ing to the bend-frequency, scales approximatel y linearly 
with black hole mass fr om AGN to GB Hs (McHardv 
19881: lEdelson fc Nandral Il999l: lUttlev et all |2002| . I2OO5F 
Markowitz et all 120031 ; iMcHardv et all |2004 120051 ). albeit 
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with some scatter. However, the scatter is entirely accounted 
for by variations in accretion rate, allowing scaling between 
AGN and GBHs on time-scales from ~ years to ~ ms 
l|McHardv et alj|2006l h 

GBHs are observed in a number of distinct X-ray spec- 
tral states which also have distinct X-ray timing proper- 
ties. Two common states are the low/hard (hereafter 'hard') 
and high/soft (hearafter 'soft') states. In the hard state, the 
energy-spectrum is dominated by a highly variable power 
law component and the PSDs are well fitted by multiple 
broad Lorentzians. For use in AGN, where signal/noise is 
lower than in GBHs, this PSD shape can be approximated by 
a double-bend power law with slopes a = 0, -1 and -2, from 
low to high frequency, where the high- and low-frequency 
bends correspond to the strongest peaks in the Lorentzian 
parameterisation. The breaks are typically separated by only 
one to two decades in frequency. In the soft state, the en- 
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ergy spectrum is dominated by an approximately constant 
thermal disc component which extends into the X-ray band 
in GBHs but which in AGN is shifted down to the opti- 
cal/UV band. Therefore, a meaningful comparison between 
the PSDs of soft state GBH and AGN can only be made in 
cases where the GBH power-law emission is strong enough 
to show significant variability. Such GBHs are rare and the 
best example is Cyg X-l which shows a '1 /f PSD over many 
decades of frequencies (|Reig et al.ll2002h . The soft state is 
distinguished from the hard state by having only one, high 
frequency, break in this power law, from slope -1 to -2. 

It has been suggested that this pure simple broken or 
cut-off power-law PSD shape is unique to the soft state of 
Cyg X-l, which is a persistent source. However in transient 
GBHs with similar X-ray spectra, the power law PSD com- 
ponent m ay be seen in combinatio n with broad Lor e ntzian 
features jPone fc Gierlih ski 2005). lAxelsson et all (120061 ) 
also note that a mixed power law plus Lorentzian PSD is 
also present in Cyg X-l in lower luminosity, harder spectral 
states, but as the luminosity rises the Lorentzian features 
weaken and the power law PSD component strengthens un- 
til, in the softest state, it completely dominates. Since the 
softest spectral states of transient GBHs are dominated by 
constant disc emission we cannot determine whether they 
show a similar PSD shape to Cyg X-l. 

However, a direct comparison of transient GBHs and 
Cyg X-l is complex, since the transients show much larger 
luminosity changes, and complex hyste resis effects in spec- 
tral hardness vers us luminosity (e.g. iHoman et al.l l200ll, 
iBelloni et ai1l2005l ) which are not seen in Cyg X-l. There- 
fore it is not clear that one can compare timing properties 
between Cyg X-l and transient GBHs simply as a function 
of observed X-ray spectrum. 

None the less, it is still interesting that the X-ray spec- 
trum of Cyg X-l never becomes totally disc-dominated, 
and always contains a relatively strong variable component 
whose PSD resembles that of X-ray bright AGN. If variabil- 
ity originates, at least partly, in the disc, so power spectral 
shape is related to the disc structure, that structure might 
be severely disrupted during outbursts, thereby suggesting 
a possible difference between the persistent Cyg X-l and the 
transient GBH sources. The similarities between the PSDs 
of Cyg X-l and AGN may also be related to the possible 
similariti es in accretion flows betw een AGN and Cyg X-l 
noted bv lDone fc Gierlinskil (|2005h . 

To date, NGC 3783 and the Narrow Line Seyfert 1 
Galaxy (NLS1) Ark 564 are the only AGN with suggested 
second, low-frequency breaks in their PSDs (i.e. similar to 
low/hard GBHs) and are both commonly referred to as be- 
ing unusual (e.g. iDone fc Gier liriski 2005|). The power spec- 
tral evidence f or a second break in the case of Ark 564 
is very strong dPounds et all l200ll ; iPapadakis et al.1 2002; 
iMarkowitz et al.ll2003l . M c Hardv et al. in prep.). Of all the 
AGN with good timing data, Ark 564 shows the highest 
accretion-rate (possibly super- Eddington) so it would not be 
surprising if it were in an unusual state, e.g. the 'very high' 
state where the PSD, in GBHs, also displays two distinct 
breaks. The properties of NGC 3783, on the other hand, are 
similar to those of AG N with proven soft-st ate PSDs (e.g. 
NGC 3227, NGC 405 1 IMcHardv et"ai]|2004 MCG-6-30-15 
IMcHardv et~al| 120051) . and in particular it is radio quiet 
(e.g. Revnoldsf 19971 ). In the hard state, GBHs are strong 



radio sour ces whereas in the soft state the r a dio emission is 
quen ched l|Corbel et all 120001 ; iFenderl l200ll ; iKording et ail 
l2006h . We also note that NGC 3783 has a more moder- 
ate accretion rate than Ark 564 (~ 7%), and more similar 
to the other AGN mentioned above, and Cyg X-l changes 
from the hard to the soft state at a round 2% of the Edding- 
ton accretion rate (i.e. itie =0.02) ( Pottschmidt et al"1l2003l ; 
I Wilms et aill2006l ; lAxelsson et al.ll2006l ). These two facts do 
not lie easily with a hard state identification of NGC 3783. 
Thus it would be surprising, and might indicate that our cur- 
rent ideas regarding the scaling between AGN and GBHs are 
not entirely correct, if NGC 3783 were proven to have a hard 
state PSD. It is therefore important to determine whether 
NGC 3783 does have a second, low frequency, break in its 

PS D or not. 

IMarkowitz et al.1 l|2003l ) recognised the presence of a 
break in the 2-10 keV PSD of NGC 3783 at 4 x 10~ 6 Hz 
and found provisional evidence for a second lower-frequency 
break at ~ 2xl0 -7 Hz. Specifically, Ma rkowitz et al.1 (2003) 
rejected the possibility that the PSD is described by a single- 
break power law with low-frequency slope -1, similar to other 
AGN, at the 98% confidence level. In this paper we re- 
investigate the evidence for the second break in the PSD of 
NGC 3783, using new long-term monitoring data that cov- 
ers the frequency range where the break appears to be. By 
including additional RXTE archival data spanning several 
years, along with short time-scale observations by XMM- 
Newton, we will demonstrate that the improved PSD is 
perfectly compatible with a single-bend power law, consis- 
tent with the behaviour of the other moderately-accreting 
Seyferts. In Section 2 we describe the observations and the 
methods by which we extract the RXTE and XMM-Newton 
light curves. In Sections 3 we discuss the PSD of NGC 3783 
as produced from the RXTE and XMM-Newton observa- 
tions, and compare it with various PSD models. In Section 
4 we briefly review the implications of our observations. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 RXTE Data Reduction 

From 1999 to 2006, NGC 3783 has been the target of various 
monitoring campaigns with RXTE. These campaigns have 
consisted of short, ~1 ks duration , observations wit h the 
proportional counter array (PCA. [Zhang et al.lll993l ). We 
have analysed the archival PCA STANDARD-2 data and 
our own proprietary data with FTOOLS v6.0.2 using stan- 
dard extraction methods. We use data from the top layers of 
PCUs and 2 up to 2000 May 12 and only top layer PCU 2 
data from observations after this date. The remaining PCUs 
were not used due to repeated breakdowns. 

Data were selected according to the standard 'good- 
time' criteria, i.e. target elevation < 10°, offset pointing 
< 0.02°, and electron contamination < 0.1. The background 
was simulated with the L7 model for faint sources using 
PCABACKEST v3.0. The response matrices for each PCA 
observation were calculated using PCARSP vlO.l. The fi- 
nal 2-10 keV fluxes were calculated using XPSEC vl2.2.1 
by fitting a power law with galactic absorption to the PHA 
data. 

The RXTE data used in our analysis, together with the 
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Figure 1. RXTE long-term light curve of NGC 3783 in the 2- 
lOkeV band. 
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Figure 3. XMM-Newton light curve in the 0.2-2 keV band of 
NGC 3783. 
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Figure 2. RXTE intense sampling light curve in the 2-10 keV 
band of NGC 3783. 



sampling patterns, are listed in Table [T] and displayed in 
Fig. [1] The early data (to MJD 52375) with 4 d sampling, 
together with the 20 d period of 3 h sampling already pre- 
sented bv lMarkowitz et all (|2003T ) are followed, after a 2 year 
gap, by our new long-term monitoring, with 2 d sampling. 
As the gap is large compared to the duration of each moni- 
toring campaign we will include data from each monitoring 
campaign as separate lightcurves in our fits. 



2.2 XMM-Newton Data Reduction 

NGC 3783 was observed by XMM-Newton during revolu- 
tions 371 and 372, between 2001 December 17 and 2001 De- 
cember 21 . Temporal analys is of these data were first pre- 
sented by iMarkowitd (|2005T ) who discusses the coherence, 
frequency-dependent phase lags, and variation of high fre- 
quency PSD slope with energy. Here we use these data to 
constrain the high frequency part of the overall long and 
short timescale PSD. We used data from the European Pho- 
ton Imaging Cameras (EPIC) PN and MOS2 instruments, 
which were operated in imaging mode. MOS1 was operated 
in Fast Uncompressed Mode and we do not use those data 
here. The PN camera was operated in Small Window mode, 
using the medium filter. Source photons were extracted from 
a circular region of 40" radius and the background was se- 



lected from a source-free region of equal area on the same 
chip. We selected single and double events, with quality 
flag=0. The MOS2 camera was operated in the Full Win- 
dow mode, using the medium filter. We extracted source 
and background photons using the same procedure as for 
the PN data and selected single, double, triple and quadru- 
ple events. These data showed no serious pile-up when tested 
with the XMM-SAS task epatplot. 

We constructed light curves, for each detector and or- 
bit, in the 0.2-2, 2-10 and 4-10 keV energy bands. We filled 
in the ~ 5 ks gap in the middle of orbit 371 light curves, and 
some other much smaller gaps, by interpolation and added 
Poisson noise. The resulting PN and MOS2 continuous light 
curves were then combined to produce the final light curves 
for each orbit. The combined, background subtracted, av- 
erage count rates in the 0.2-10 keV band were 11.8 c/s for 
orbit 371 and 15.8 c/s for orbit 372, and the 0.2-2 keV com- 
bined light curve is shown in Fig. [3] Poisson noise dominates 
the PSD on timescales shorter than 1000s, so the light curves 
were binned into 200s bins. 



3 POWER SPECTRAL MODELLING 

3.1 Combining RXTE and XMM-Newton data 

To determine the PSD over the largest possible frequency 
range we combine the RXTE and XMM-Newton data. In 
GBHs the break-frequency and slope of the PSD below the 
break appear to be independent of the chosen energy band 
ilCui et all [l997l: JChurazov et all l200ll; iNowak et all Il999l : 



iRevnivtsev et al-lboOOl ; iMcHardv et al]|2004h . On the other 
hand, the PSD normalisation and the slope abo ve the break 
are often energy-dependent ( Markowit3 I2005T ). Therefore, 
when combining data from different instruments, it is prefer- 
able to use similar energy ranges. The RXTE data are in 
the 2-10 keV band and, for NGC 3783, that band has a 
median photon energy of 5.7 keV. The XMM-Newton band 
with the same median photon energy is 4.1-10 keV. How- 
ever the count rate in that XMM-Newton band is low (2 
c/s) so we only detect significant source power above the 
Poisson noise level at frequencies below 10 _4 Hz. To probe 
higher frequencies we can use the 0.2-2 keV XMM-Newton 
data (8.8 c/s) but we must re-scale its PSD normalisation to 
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Light curve 


Sampling interval 


Observation length 


Date Range [MJD] 


RXTE Long-term 1 


~4.36 days 


1194.6 days 


51180.5-52375.1 


RXTE Long-term 2 


~2.1 days 


928.3 days 


53063.4-53991.6 


RXTE Intense monitoring 


~3.2 hours 


19.9 days 


51960.1-51980.1 


XMM-JVewton observations (2 orbits) 


200-s 


3.2 days 


52260.8-52264.0 



Table 1. Summary of the RXTE and XMM-JVewton light curves used in the analysis of NGC 3783, including their sampling frequency 
and date range. 



that of the 4-10 keV PSD. We determined the scaling cor- 
rection by producing PSDs in both energy bands and fitting 
the same bending power law model to the noise-subtracted 
data. On the assumption that the PSD shape below the high 
frequency break is energy-independent, the combined RXTE 
2-10 keV and XMM-Newton 0.2-2 keV PSD will then have 
the shape of the 0.2-2 keV PSD. 



3.2 Monte Carlo simulations 

We use the Monte Carlo technique of lUttlev et all (|2002T I 
(PSRESP), to estimate the underlying PSD parameters in 
the presence of sampling biases. In this method we first 
calculate the observed (or 'dirty') PSD, in parts, from the 
observed lightcurves, using the Discrete Fourier Transform. 
Here the PSD estimates are binned in bins of width 1.3u, 
where v is the st arting frequency by taking th e average of 
the log of power ()Papadakis fc Lawrencdll993l L We require 
a minimum of 4 PSD estimates per bin. We then compare si- 
multaneously the dirty PSDs from each lightcurve with var- 
ious model PSDs derived from lightcurves simulated with 
the same sampling pattern as the real observations. We al- 
ter the model parameters to o btain the best fit for any given 
model. We refer the reader to lUttlev et all <|2002h for a full 
discussion of the method. 

For each set of chosen underlying-PSD model param- 
eters, we simulate red-n oise light curves, as described by 
iTimmer fc Koenid (Il995l ). The RXTE light curves are sim- 
ulated with time resolutions of 10.5 h, 5.0 h and 18 m for the 
first and second long time-scale and the medium time-scale 
light curves respectively. The simulated resolution, which is 
10 times shorter than the typical sampling intervals of the 
real observations, given in column 2 of Table[T] is to take into 
account the effect of aliasing. These simulated light curves 
were resampled and binned to match the real NGC 3783 ob- 
servations. XMM-Newton light curves were simulated with 
200-s resolution, as at shorter time-scales the underlying 
varability power is negligible compared to the Poisson noise, 
so aliasing does not play a role. The Poisson noise level was 
not subtracted from the observed PSD but was added to 
the simulated PSDs. To reproduce the effect of red-noise 
leak, each light curve was simulated to be ~300 times longer 
than the real observation, and was then split into sections, 
constituting 300 simulated light curves for each observed 
lightcurve. The simulated model average PSD is evaluated 
from this ensemble of PSD realisations, and the errors are 
assigned from the rms spread of the realisations within a 
frequency bin. 

We present the results of several PSD model fits in an 
attempt to quantify the underlying model shape that best 
describes the PSD of NGC 3783, and associate an aCCep- 
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Figure 4. The best fit unbroken power law PSD of the combined 
RXTE and XMM-Newton data. The solid lines represent the ob- 
served data and the points with error bars exhibit the biased 
model and the spread in individual realisations of the model. The 
dashed line is the underlying model used to generate the simu- 
lated PSD. The three lowest frequency data sets are from RXTE 
observations and the high frequency data set (~10 — 5 Hz) is from 
the XMM-Newton observations. Note that the rise in power at 
the highest frequencies is due to the photon Poisson noise. 

tance probability with each model. We initially test a sim- 
ple unbroken power law model. We next fit a power law with 
a single-bend in the PSD, and then a model incorporating 
a double-bend. We also fit a single-bend power law with a 
Lorentzian component. 

3.3 Unbroken power law model 

To begin, we fitted a simple power law model to the data of 
the form: 

where A is the normalisation at a frequency vq, and a 
is the power law slope. We made 900 simulations and in Fig. 
[4] we show the best fit plotted in v x P v , which has a PSD 
slope of -2.1. However, the fit is poor and this model can be 
rejected with a probability > 99 %, or > 3 a. 

3.4 Single-bend power law model 

Here we fit a single-bend power law to the data. This model 
best describes the PSD of Cyg X-l in the high/soft state, 
and provides a goo d fit to the PSDs of the AGN NGC 4051 
and MCG-6-30-15 l|McHardv et al.ll2004l2005l ). 
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Figure 5. The best fit single-bend power law PSD of the com- 
bined RXTE and XMM-Newton data. The various lines represent 
the same data as seen in Fig. [4] 



t — i — i — i — I 




\ 



| Q I — — — — — i i — 

-10~ 8 10" 7 10" 6 10" s 10" 4 10" 3 

Frequency (Hz) 

Figure 6. The best fit double-bend power law PSD of the com- 
bined RXTE and XMM-JVewton data. The various lines represent 
the same data as seen in Fig. [J] 
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Figure 7. Single-bend power law model: 68, 90, and 99% confi- 
dence contours for the bend frequency, i/g , and the high frequency 
slope, —an, for the single-bend power law fit to the combined 
RXTE and XMM-Newton PSD. 
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Figure 8. Single-bend power law model: 68, 90, and 99% confi- 
dence contours for the bend frequency, ub, and the low frequency 
slope, — DLL, for the single-bending power law fit to the combined 
RXTE and XMM-Newton PSD. 
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Fig. [5] presents the observed PSD fitted with a single- 
bend power law model, for which a good likelihood of accep- 
tance is obtained (P = 44 %). The best fit bend- frequency is 
vb = 6.2lg°g 6 x f0~ 6 Hz, the high-frequency slope is an = 
— 2.6l°' 6 , and the low-frequency slope is ai = — 0.8^0 g- 
The errors are 90% confidence limits, an asterisk indicates 
that the limit is unconstrained. For oih the best fit value 
is well within the searched parameter space but the degen- 
eracy produced by red-noise leak in the probability at high 
values of an, means that the upper limit is not constrained 
at the 90% confidence level. The confidence contours for the 
main interesting parameters are plotted in Figs. [7] and [S] Ta- 
ble[2]shows the single-bend power law best fit parameters to 
the data. The best fit single-bend f requency obtained here i s 
consistent with the value found by iMarkowitz et al.l (|200of ) 

3.5 Double-bend power law model 

IMarkowitz et all (|2003h provide tentative evidence that 



a second, lower, frequency break exists in the PSD of 
NGC 3783. Thus, we also fitted a more complex double-bend 
power law model to see if the goodness-of-fit is improved sig- 
nificantly. The double-bend power law model is given by: 



Av° 



P{y) = 



where oti is the intermediate-frequency slope and vl 
and vh are the low and high bend- frequencies respectively. 
We fixed the low-frequency slope to zero, to avoid making 
the simulations computationally prohibitive, and because a 
low-frequency slope of zero would allow the best qualitativ e 
comparison to the low state of Cyg X-l l|Nowak et al.|[l999l ). 

Fig. [6] presents the same observed PSD as in Fig. [5] 
but fitted with the double-bend power law model. A good 
likelihood of acceptance is obtained (P=64 %). The best- 
fitting high bend-frequency is vu = 2.6^* x 10~ 5 Hz, the 
high-frequency slope is an = — 3.2^J;' 2 , the intermediate- 
frequency slope is qj = — 1.3l*, the low-frequency bend is 
vl = 1-71* x 10 -7 Hz. As before, we use 90% confidence lim- 
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Model 


Normalisation 

(a) 




a i 


OIL 


''// 
(Hz) 


(Hz) 


Acceptance 

(96) 


Single-bend 
Double-bend 


1.5 x icr 4 

1.0 x 10 2 


-2 6 +0 - 6 
-3.2±i' 2 


NA 

-i.3i; 


n s+ - 8 
0.0 


6.2±f 6 6 x lO" 6 
2.6±* x 10~ 5 


NA 
1-Tt* x 


44.4 
63.9 



Table 2. Best fit model parameters for the examined models to the combined RXTE and XMM-Newton PSD of NGC 3783. The errors 
on the single- and double-bend fits are calculated from the 90% confidence intervals. The bend-frequency for the single-bend model, ub, 
is denoted here as vjj . An asterisk indicates that the limit is unconstrained. 



its. The added parameters allow extra freedom to find better 
fit probabilities for any given set of double-bend parameters. 
For this reason, the contour levels cover larger ranges in the 
parameter space and therefore, most of the 90% contours 
in our double-bend fit remain unbounded over the fitted 
parameter space. The high-frequency slope is subject to the 
same problems as in the single-bend model. Table[2]contains 
a summary of the best-fitting model parameters. 

The best-fitting low-frequency bend is found close to the 
lowest frequencies probed by the data and, as seen in Fig. [9] 
it is essentially unbounded down to the lowest measurable 
frequency at the 68% confidence level. These facts suggest 
that the second, low-frequency, bend might not be required 
by the data and that the improvement in the fit might be 
only due to the increased complexity of the model fitted. 

The likelihood of acceptance is better in the double- 
bend model than in the single-bend model, 64 versus 44 % 
respectively, but there are more free parameters. In order 
to determine the significance of this improvement, we per- 
formed the following test. Using the best-fitting single-bend 
PSD parameters, we generated 300 realisations of the sets of 
RXTE and XMM-Newton lightcurves. Each realisation was 
then fitted with the best-fitting double-bend parameters, ex- 
actly as was done with the real data, and the distribution 
of their fit probabilities was constructed. We found that 121 
out of the 300 single-bend simulations have a higher fit prob- 
ability than the real data, when fitted with the double-bend 
model. Therefore, we conclude that the improvement in fit 
probability is no more than may be expected from fitting a 
model which is more complicated than required by the data: 
the double-bend model does not represent a significant im- 
provement. 



3.6 Single-bend power law with a Lorentzian 
component 

We finally consider whether the observed PSD might be 
best-described by adding a Lorentzian component, such as 
are commonly used to describe broad-band noise compo- 
nents in GBHs (e.g. lNowakll2000h . to the single-bend power 
law. We are motivated to consider this possibility because 
the PSD of the intense-sampling RXTE light curve is not 
very well described by either the single- or double-bend 
power law model. Visual inspection of this light curve, shown 
in Fig. [2] suggests that the variability is strongly concen- 
trated on time-scales of around a day, or equivalently, fre- 
quencies around 10 Hz, which is confirmed by the peak 
seen in the corresponding section of the PSD, and the drop 
in the same PSD at lower frequencies (~ 10 _6 Hz). The 
long-term monitoring PSDs, however, do not show a dip 
at 10 _6 Hz, creating a large discrepancy in the PSD mea- 




Figure 9. Double-bend power law model: 68, and 90% confi- 
dence contours for the high bend-frequency, vh , and the low bend- 
frequency, vl, for the double-bend power law fit to the combined 
RXTE and XMM-Newton PSD. 



surements at this frequency. A strongly peaked component 
in the underlying PSD, at ~ 10~ 5 Hz, could produce the 
observed features. Such a component would appear as a 
peak in a PSD that covered frequencies above and below 
its peak-frequency, but would be insufficiently sampled by 
the long-term monitoring campaigns; thus, its power would 
be aliased into the highest frequencies of the longer time- 
scale data, making them rise above the underlying model 
level and causing the apparent disparity. 
The Lorentzian profile is described by: 



AQv c 



+ 4Q 2 (iv 



where the centroid frequency u c is related to the peak- 
frequency v p by v p — v c yT+ 1/4Q 2 and the quality factor 
Q is equal to v c divided by the full width at half maximum 
of the Lorentzian. The variable A parameterizes the relative 
contribution of the power law and Lorentzian components 
to the total rms. Fitting a Lorentzian component in addition 
to a single-bend power law provides a good fit (P=52 %). 
The best-fitting Lorentzian contributes 20% of the variance 
in the frequency range probed and its best-fitting parame- 
ters are quoted in Table [3] Fig. HOI shows the observed PSD 
compared with the best-fitting single-bend power law model 
plus a Lorentzian component. The Lorentzian feature in the 
model can reproduce qualitatively the spurious power at the 
high frequency end of the long-term monitoring data and the 
turn down effect observed in the intensive-sampling data. 

To determine the significance of the Lorentzian compo- 
nent fit we repeated the procedure used in determining the 
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(Hz) 


(Hz) 


Q 


A 




Off 


Acceptance 

(96) 


4-8t5.g X iO" 6 


1.1±°;| x 10" 5 


°- 1 -3.6 


0.9± ( * 7 


-l-0±* 


-2.6±* 


52.3 



Table 3. Best fit single-bend power law with Lorentzian component model parameters to the combined RXTE and XMM-Newton PSD 
of NGC 3783, where v v is the Lorentzian peak frequency, Q is its quality factor, ug is the power law bend frequency and oil and ajj 
are the power law slopes below and above the bend, respectively. The errors are calculated from the 68% confidence intervals, and an 
asterisk indicates that the limit is unconstrained. 
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Frequency (Hz) 

Figure 10. The best-fitting single-bend power law with a 
Lorentzian component. The fit was done using the entire data 
set but here we only show the Lorentzian region. As before, solid 
lines represent the real data PSD, dashed lines represent the best- 
fitting model and markers with error bars represent the model 
distorted by sampling effets. The Lorentzian feature in the model 
can reproduce qualitatively the spurious power at the high fre- 
quency end of the long-term monitoring data and the turn down 
effect observed in the intensive-sampling data. 



significance of the double-bend model. We found that 222 of 
the 300 single-bend simulated PSDs have a higher fit proba- 
bility than the data, when fitted with the single-bend power 
law plus Lorentzian model. This result indicates that the 
increase in fit probability could be due to the added com- 
plexity of the model, and that the improvement in the fit 
over a simple bending power law is not significant. 



4 DISCUSSION AND CONCLUSIONS 

We have combined our own new RXTE monitoring data 
with archival RXTE and XMM-Newton observations to 
construct a high-quality PSD of NGC 3783 spanning five 
decades in frequency. 

We find that a 'soft' state model, with a single 
bend at 6.2 x 10~ 6 H z, similar to that found earlier by 
iMarkowitz et all (120031 ). a power law of slope approximately 
-0.8 extending over almost three decades in frequency below 
the bend, and slope above the bend of approximately -2.6 
is a good fit to the data. We also find that a 'hard' state 
model, with a double bend, fits the data, as does a model 
with a single bend plus an additional Lorentzian compo- 
nent. However the improvement in fit is marginal and, given 
the additional free parameters, is not significant. Thus we 
conclude that a simple 'soft' state model provides the most 
likely explanation of the data. 



Assuming a ma ss of 3 x 10 7 A/ o for NGC 3783 
jPeterson et al.l l2004h , and an accretion rate of 7% of 
the Eddington lim it ijllttlev fc McHardvl 120051 . based on 

1 Woo fe Urrvll2002l ), then NGC 3783 is still in good agree- 
ment with the scaling of PSD brea k timescale as ~ M/rh g 
between AGN and GBHs found bv lMcHardv et all (|2006D . 

Our new fits, show that the PSD of NGC 3783 is per- 
fectly consistent with a single-bend power law with low- 
frequency slope o f -1, i n contrast with the earlier result of 
IMarkowitz et all (|2003T l. who found that a similar model 
was rejected tentatively at ~ 98% confidence. The differ- 
ence can be understood in terms of the improved long-term 
data. Our new RXTE monitoring observations occur every 2 
days, compared to 4 days previously, thereby increasing the 
long term RXTE data set by a factor 2.6 and, in particular, 
providing overlap at high frequencies with the RXTE inten- 
sive monitoring data. The drop in long-timescale variability 
power, evident in the older long term monitoring data is not 
reproduced by the new long-term monitoring data, showing 
that this drop could be just a statistical fluctuation. In addi- 
tion, the very high frequencies are better constrained by the 

2 orbits of X MM-Newton data than by the earlier Chandra 
data used by IMarkowitz et alj (|2003l ). 

The classification of the PSD as being 'soft' state means 
that NGC 3783 is no longer considered unusual amongst 
AGN. The fact that this AGN is radio-quiet strongly sup- 
ports the analogy with GBHs in the so ft state. Also the ac- 
cretio n rate of N GC 3783 ( m E =0.07 ) l|Uttlev fe McHardvl 
120051 . based on IWoo fc Urrvl I2002T I is similar to that 
of other AGN with soft-state PSDs (e.g. NGC 3227 



lUttlev fe McHardvl 120051. NGC 405 1 iMcHardv et"ai] |2004| . 
MCG-6-30-15 iMcHardv et al.ll2005l ). This accretion rate is 
above the rate at which the persistent GBH Cyg X-l tran- 
sits between hard and soft states in either direction and at 
which o ther GBHs transit from the soft to hard state ( vtie 
=0.02) (|Maccarone et al.ll2003l ; lMaccaronell2003l ). We note 
that other transient GBHs in outburst, where the variable 
power law emission in the soft state PSD is weak, can re- 
main _^n_tli£l2^rd_st^te tojnuch higher accretion rates (~ 2- 
50% 1 oman fe Bellonill2005l ) but it is not clear whether we 
should expect similar PSD shapes to AGN for such outburst- 
ing sources. Thus NGC 3783 remains compatible with other 
moderately accreting AGN in being analogous to Cyg X-l 
in the soft state. It is, of course, possible that the transition 
rate might not be independent of mass. Observations do not 
yet greatly constrain the transition rate as a function of mass 
but the abscence of large deviations from the so-called 'fun- 
damental' plane of radio luminosity , X-ray luminosity and 
black hole mass l|Merloni et al.ll2003l ; iFalcke et ai1l2004l ) ar- 
gues aga inst a large spread in the transistion accretion-rates 
(e.g. see lKording et al]|2006l ). In the case of Seyfert galaxy 
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NGC 3227, the a ccretion rate is ~ 1-2% a nd a 'soft' state 
PSD is measured l|Uttlev fc McHardvl[2005h . which suggests 
that the transition accretion-rate in AGN should be at or 
below that value. 

Our observations, which show that NGC 3783 does not 
have a highly unusual PSD, therefore confirm the growing 
similarities between AGN and Galactic black hole systems 
and leave only Arakelian 564, which is probably a very high 
state object, as the only AGN showing clear double breaks 
(or m ultiple Lorentzians) in its PSD (e.g. lArevalo et al.l 
120061 . M c Hardy et al. in prep.). 
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